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An important objective of the present joint project with Atmospheric and Environ- 
mental Research, Inc., is the development of a couple 2-D model of dynamics, radiation 
and chemistry. Our existing model (see Ko et al. (1985)) is not fully coupled, in the sense 
that net radiative heating, Q, and the isentropic eddy diffusion coefficient, K vv , have to 
be separately specified. 

Recent developments (Plumb and Mahlman, 1986; NASA/WMO, 1985) suggest that in 
the geostrophic limit, the two quantities may be related. It is, however, rather problematic 
to implement such a geostrophic result in a zonally averaged model because (i) the model 
is global, while geostrophy breaks down near the equator, and (ii) geostrophy filters out 
equatorial waves and gravity waves, which are known to be important components in 
the eddy forcing of the mean flow. Furthermore, the common procedure of adopting the 
geostrophic result is based on the assumption of conservation of the quasi-geostrophic 
potential vorticity along isobaric surfaces. (Newman et al., 1986). It is known that a 
better conserved quantity is the (Ertel’s) isentropic potential vorticity along isentropic 
surfaces. 

Our effort during the first six months of the project is to formulate a nongeostrophic 
theory that is more compatible with our global 2-D model. By relating the Eliassen-Palm 
flux divergence, which represents the net eddy forcing of the mean zonal momentum, to 
the flux of the isentropic potential vorticity, we have been able to deduce a nongeostrophic 
relationship between Q and K vv in isentropic coordinates. This work is documented in 
Tung (1986), submitted for publication in J. Atmos. Sci.. 

We have also been able to infer seasonal and latitudinal distributions of K yv from the 
mean momentum budget. The result will be submitted for publication shortly. By adopt- 
ing a nongeostrophic theory and using Ertel’s isentropic potential vorticity in isentropic 
coordinates instead of quasi-geostrophic potential vorticity in pressure coordinates, many 
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of the large-scale negative K vv regions encountered by Newman et al. (1986) disappeared. 
The remaining large region of negative K vy occurs in the easterly region of the summer 
hemisphere, but the magnitudes of K vv there are extremely small due to the inability 
of the planetary-stationary waves to significantly penetrate into this easterly region in 
stratosphere. 

During the next six months, an algorithm based on the above-mentioned nongeostrophic 
relationship will be implemented into the AER 2-D model. This will remove K yv as an 
independent input parameter for the model. Self-consistent calculations will be performed 
for various chemical species, including ozone. 
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Nongeostropic Theory of Zonally Averaged 
Circulation 
Part I: Formulation 


Ka Kit Tung 

Department of Mathematics, 
M.I.T., Cambridge, MA 02139 


ABSTRACT 


A nongeostrophic theory of zonally averaged circulation is formulated 
using the" nonlinear primitive equations on a sphere, taking advantage of 
the more direct relationship between the mean meridional circulation and 
diabatic heating rate which is available in isentropic coordinates. Possible 
differences between results of nongeostrophic theory and the commonly used 
geostrophic formulation (e.g. Edmon et al (1980)) are discussed concerning 
(a) the role of eddy forcing of the diabatic circulation, and (b) the “non- 
linear nearly inviscid” limit (Held and Hou, 1980) vs the geostrophic limit. 
Problems associated with the traditional Rossby number scaling in quasi- 
geostrophic formulations are pointed out and an alternate, more general 
scaling based on the smallness of mean meridional to zonal velocities for a 
rotating planet is suggested. Such a scaling recovers the geostrophic balanced 
wind relationship for the mean zonal flow but reveals that the mean merid- 
ional velocity is in general ageostrophic, analogous to the semi-geostropic 
theory of Hoskins and Bretherton (1972). 

A set of general diagnostic tools comparable in scope to their geostrophic 
counterparts is given in Part I, including (a) a generalized definition of 
Eliassen-Palm flux divergence (without restriction to small amplitudes, steady 
state or to adiabatic flows), the vanishing of which is a necessary condition 
for nonacceleration (Andrews, 1983), (b) a nonlinear formula that relates the 
flux of Ertel’s potential vorticity to the Eliassen-Palm flux divergence and 
(c) a relationship between the Eliassen-Palm flux divergence and isentropic 
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mixing coefficient, K vv , used in chemical tracer transport equations in isen- 
tropic coordinates. From the mean momentum budget, we give in Part II an 
estimate of the Eliassen-Palm flux divergence using fitted “observed” field 
of net radiative heating rate. From this an estimate of the magnitude and 
latitudinal/seasonal variation of K vv is also provided. 
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1. Introduction 


Great progress has been made in our understanding of the zonally aver- 
aged circulation in the atmosphere since the introduction of the concept of 
residual circulation by Andrews and McIntyre (1976) and Boyd (1976), which 
helps put in a more proper perspective the role of “eddies” (i.e. deviations 
from zonal symmetry) in driving the mean (zonally symmetric) circulation. 
However, our present understanding is based largely on a geostrophic ver- 
sion of the general theory (see Edmon et al (1980), Dunkerton et al (1981), 
Palmer (I981ab) and Andrews et al (1983)). In its geostrophic form, the set 
of transformed zonally averaged equations of motion shows clearly that, in 
the absence of eddies in the form of an Eliassen-Palm flux divergence, eddies 
do not accelerate the zonal mean flow, a result first recognized by Charney 
and Drazin (1961) and Dickinson (1969). A more controversal consequence of 
the geostrophic theory is that at equilibrium, an “almost frictionless” (molec- 
ular diffusion only say) stratosphere in the absence of large-scale eddies would 
be “extremely close” to radiative equilibrium, in which the absorption of solar 
insolation is simply balanced by an increase in local temperature, without 
inducing a meridional circulation (see Mahlman et al (1984)). It is therefore 
often concluded (see WMO/NASA (1985), Chapter 7) that global pattern of 
rising and descending motions in the stratosphere “owes its existence to the 
presence of asymmetric motions”. 

The above cited conclusion appears to have been contradicted by a num- 
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ber of zonally symmetric calculations based on the nonlinear primitive (i.e. 
nongeostrophic) equations that produce realistic looking Hadley circulations 
in the absence of large-scale eddies but in the presence of small-scale mixing 
(i.e. viscosity) (Schneider and Lindzen, 1977; Schneider, 1977; Nakamura, 
1978; Held and Hou, 1980). We can further show that even the small-scale 
eddies are not needed by demonstrating that global scale nonlinear out-of- 
radiative equilibrium circulations exist in both inviscid and “nearly inviscid” 
limits, provided that nongeostrophic equations are used. 

While these “thought experiments” of a hypothetical axisymmetric at- 
mosphere are useful in highlighting the qualitative difference between the 
geostrophic and nongeostrophic formulations, they do not address the more 
practical question: Is our atmosphere in the geostrophic regime? Or is the 
“nonlinear nearly inviscid” regime (Held and Hou 1980) more applicable? 
These issues cannot be addressed using a formulation that adopts the a pri- 
ori assumption of geostrophy. 

It is important to recognize that there is no a priori justification for the 
geostrophic scaling when applied to the zonally averaged zonal momentum 
equation. When eddy forcing is absent, the geostrophic form of the mean 
zonal momentum equation suggests that the mean meridional velocity, and 
so the Coriolis acceleration, vanishes, while geostrophic scaling requires that 
Coriolis acceleration dominates relative acceleration terms. Therefore, it is 
obvious that one cannot justifiably use the geostrophic equations to deduce 
results concerning what happens in the absence of large-scale eddies. On 
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the other hand, the above observation does not necessarily imply that the 
geostrophic approximation cannot be applied to the stratosphere; given suf- 
ficiently strong magnitudes of eddy Eliassen-Palm flux divergence of global 
length scales, a scaling based on the smallness of the Rossby number may in- 
deed be valid. It is therefore an objective of the present work to (i) point out 
the different regimes an atmosphere can be in depending on the magnitudes 
of eddy forcing, and (ii) to find out which of these regimes our atmosphere is 
in during different seasons. Task (i) is discussed in the present Part I, while 
task (ii) is relegated to Part II. 

Another objective of the present work is to develop a complete set of 
general diagnostic tools that are comparable in scope to the correspond- 
ing geostrophic diagnostics. 1 In order to make our diagnostics more easily 
adaptable to current procedures for data analysis, we have emphasized the 
derivation of relationships of various eddy diagnostic quantities to Ertel’s 
potential vorticity (Ertel, 1942). McIntyre and Palmer (1983, 1984) have ef- 
fectively argued for the usefulness of Ertel’s potential vorticity on isentropic 
surfaces as a diagnostic tool for visualizing large-scale nonlinear dynamical 
processes in the stratosphere. Hoskins et al (1985) have reviewed the use and 
significance of isentropic potential vorticity maps in a variety of atmospheric 

situations, not necessarily rest ricted to the stratosphere. In particular, they 
1 Such a set of diagnostics for nongeostrophic flow is hitherto not available, although a general 

(nongeostrophic) definition of the Eliassen-Palm flux divergence in pressure coordinates has 

been given by Andrews and McIntyre (1976, 1978a) and Boyd (1976) and used by Andrews 

et al (1983). 
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emphasized the Lagrangian conservation principle for ErtePs potential vortic- 
ity, which presumably is better conserved than that of the quasi-geostrophic 
potential vorticity along isobaric surfaces (Charney and Stern, 1962). In 
addition to the two main “principles” associated with the use of isentropic 
potential vorticity (namely, the conservation principle just mentioned, and 
the “inverterbility” principle (see Hoskins et al (1985)), we add further that 
(i) the use of isentropic maps of potential vorticity allows us to diagnos- 
tically calculate the Eliassen-Palm flux pseudo-divergence, which represents 
the net eddy forcing of the zonal mean flow in isentropic coordinates, and (ii) 
the Lagrangian conservative properties of the isentropic potential vorticity 
endows it with the property of a passive tracer; the dual (dynamical and pas- 
sive) character of the isentropic potential vorticity allows us, in principle, to 
deduce the transport characteristics of the atmosphere form its momentum 
budget. 

Although the relationship between the flux of geostrophic potential vor- 
ticity and geostrophic Eliassen-Palm flux divergence is well-known in pressure 
coordinates (see Edmon et al (1980), and references therein), the correspond- 
ing relationship between the flux of Ertel’s potential vorticity along isentropic 
surfaces and the Eliassen-Palm flux divergence does not appear to have been 
derived. 2 In other words, a set of general diagnostic tools for studying waves 

and wave-mean flow interactio ns is not available at the present time except 
2 However, a weak-amplitude limit of such a relationship appears to have been known to 

Andrews (1983), as can be inferred from his Eq. (4.1). 
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in the geostrophic limit in pressure coordinates. Partly as a consequence, 
Clough et al (1985) have to resort to a comparison of the isentropic maps 
of Ertel’s potential vorticity with the Eliassen-Palm flux divergence calcu- 
lated based on the quasi-geostrophic definition of Edmon et al (1980) for 
pressure coordinates. Although this somewhat inconsistent procedure is at 
present largely dictated by the retrieval process for the observational data, 
it would have been conceptually clearer if the relationship between the flux 
of Ertel’s potential vorticity along isentropic surfaces and the Eliassen-Palm 
flux divergence in the same coordinates had been known and a direct (albeit 
approximate) comparison made of these two quantities. This is especially 
true for diagnosing model generated data: Although the numerical data are 
generated using a primitive equation model, the diagnostics by Dunkerton 
et al (1981) are again based on the quasi-geostrophic formulation of Edmon 
et al (1980). 

A relationship between the flux of Ertel’s potential vorticity along isen- 
tropic surfaces and the Eliassen-Palm flux divergence is obtained without 
quasi-geostrophic approximation. This crucial relationship links the flux of 
a quasi-conservative wave property to the net wave forcing of the mean flow 
without restriction to small wave amplitudes, thus making it a useful diag- 
nostic tools applicable even to the “surf zones” (McIntyre and Palmer, 1983) 
in the stratosphere. 

The link between the net eddy forcing term in the mean zonal momen- 
tum budget and the ErteFs potential vorticity, which is quasi-conservative, 
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allows us to deduce from the mean momentum budget and mean isentropic 
gradient of potential vorticity a quantity, K yv , that appears ubiquitously 
in chemical tracer transport equations in isentropic coordinates as the isen- 
tropic diffusion coefficient (see Tung, 1982; 1984; Ko et al, 1985). Thus it 
appears that an assessment of the magnitude of the Eliassen-Palm flux di- 
vergence could also possibly afford us a preliminary look at the magnitude 
and latitudinal distribution of K VJI , that has been difficult to obtain by di- 
rect evaluation from Eulerian data on transient waves. It also turns out 
that the use of Ertel’s potential vorticity in isentropic coordinates to deduce 
K yy is less problem prone than a similar procedure in isobaric coordinates 
using quasi-geostrophic potential vorticity (Newman et al (1985), personal 
communication) . 


The use of isentropic coordinates turns out to be important in our present 
formulation of a nongeostrophic theory of zonally averaged circulation. Al- 
though many of our formulae may have their counterparts in pressure coordi- 
nates, the role played by the eddy Eliassen-Palm flux divergence in the forcing 
of the residual zonal mean circulation becomes more difficult to understand 
in pressure coordinates due to the presence of mean ageostrophic meridional 
circulations. The simple, direct relations between the Eliassen-Palm flux 

divergence and the geostrophic residual mean circulation in pressure coordi- 

\ 

nates, as discussed succinctly by Edmon et al (1980), are no longer available 
for nongeostrophic flows. Such a problem does not appear in isentropic co- 
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ordinates. For our purposes, the advantage of using potential temperature 


0 = T 



R/cp 


as the vertical coordinate (instead of pressure p, say) lies in the direct rela- 
tionship between the vertical “velocity”, 6 = f t 0, and the diabatic heating 
rate, Q, as given by the thermodynamics equation 



Thus there should be no mean circulation in the latitude-height plane for 
an atmosphere in radiative equilibrium. Conversely the presence of a meail 
meridional circulation, driven by, say, eddy forcings in the form of Eliassen- 
Palm flux divergence, necessarily induces a diabatic heating, leading to a non- 
radiative-equilibrium for the atmosphere. Such a direct relationship between 
radiation and dynamics is one of the many reasons why it is convenient to 
adopt the isentropic coordinates for our discussions to follow. The relevant 
equations will also be recast into log-pressure coordinates later in section 7. 
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2. The 2-D equations 


A listing of the set of primitive equations in isentropic coordinates can 
be found, for example, in Appendix A of Tung (1982). [There is however a 
typographical error in Eq. (A4) there, where should read cos<p]. The 
so-called 2-D equations are obtained by taking the zonal average of the 3-D 
equations by the operation 

h(<p,0) = ±-J* r h(\,<p,0)d\, ( 2 . 1 ) 

where A is the longitude, <p the latitude and 6 is the potential tempera- 
ture, used here as'the vertical coordinates. We shall use primes to denote 
deviations from zonal average, i.e. 

h! ~h~h ( 2 . 2 ) 


The primed quantities arise from asymmetric perturbations, and are referred 
to here as “eddies” or “waves” . In isentropic coordinates it is sometimes more 
meaningful to use density-weighted zonal averages, defined as (see Gallimore 
and Johnson (1981 a, b)), 

h = peh/p g , (2.3) 


where 



is the “density” in isentropic coordinates, i.e. it is mass per unit pseudo 
volume (with dz replaced by d$). The deviation from the density- weighted 
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average is denoted by an asterisk, i.e. 


h* = h - h (2.4) 

[Note that p 6 h * = 0, but h* ^ 0 in general.] 

The difference between h* and h! is ( p' e h')/~pe , which is quadratic in wave 
amplitudes. The distinction between the two is significant only for finite 
amplitude waves. But since we are interested in deriving relations that hold 
in the presence of finite amplitude perturbations, the distinction will be made 
in the present work. 

As in Tung (1982), we use capital letters to denote mass flow rates, thus 
U = pgu, V = p e vcos(p,W = ped. (2.5) 

With those definitions, the zonally averaged equations take the following 

I 

form. 


2.1. Equation of mass conservation. 


where y = a sin <p. 


k" + Ty V + Ts* =° 


( 2 . 6 ) 


2.2. Thermodynamics equation. 


The thermodynamics equation 

7t M = Q/T 
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can be written as 


W = (,/ r) = q/i 


(2.7) 


where q =' pQ is the diabatic heating rate per unit physical volume divided 
by c p , and T is the static stability parameter defined as 

r = Tde_ 

~ e dz 

The derivation of (2.7) was previously given by Tung (1982). 

[To circumvent any reference to the height coordinate z in our present formulation 
in ^-coordinate, we note that 

, 9 dz 9 9 d 

— A 

~ T dd ~ pT P6 ~ gpT d8 P ’ 

where the definition p$ = p~ and the hydrostatic relationship 

d 


have been used A further use of the ideal gas law; p = pRT, then yields 


r -l 


/fZN dlnp c p f dtnT\ 

\g ) dlnd g \ dln8 ) 


This last expression is the one we will use, assuming T = However, in 

dealing with observational data, we usually have 9 = 8(\,<p,p), in which case V can 
be obtained from 

r _ 9 dln.9 _ g_ + dT 

Rdtnp c p Hdln (m.y 

where 

H = RT/g. ) 
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2.3. Zonal momentum equation: 


d d d * 

— {p 6 u cos <p) -j- — (v'ucos^) + — (W tt cos <p) — fV = V • F, (2.8o) 

where / = 2fisin<p is the Coriolis parameter. The derivation of Eq. (2.8) is 
given in Appendix A. 

a. Eliassen-Palm flux divergence: 


Of importance in our present discussion is the right-hand side of Eq. 
(2.8a), which represents a generalized version of the so-called Eliassen-Palm 
flux divergence. Written in component form (recalling y = a simp), it is 


V • F = 


h Fy+ 


dd' 


with 

F v = —pfU*v* cos 2 <p 


and 


F $ = — cos ip 


1 d 

Ptu'O' - 

g ox 


(2.9a) 


is defined as aco f^ dx > and $ denotes the Montgomery streamfunction). 


A physical interpretation of V • F can be given, following Andrews (1983) 
but generalizing to include diabatic eddies: V • F can be viewed as the 
x-component of the force exerted by eddies on a thin tube with its axis 
oriented along the x-direction bounded by the fixed lateral sides located at 
y and y + dy, and undulating bottom and top isentropes 0 and 0 + d0. The 
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horizontal component of the Eliassen-Palm flux, F v , represents lateral flux 
of zonal momentum into the tube.' The first term in the vertical component 
Fg represents the vertical, i.e. cross isentropic, flux of zonal momentum into 
the tube, while the second term in Fg represents the x-projection of the net 
pressure force pushing against the (slanted) isentropes (the so-called form 
drag). 

The use of ( )* instead of ( )' is necessary in (2.9a) in order for (2.8a) 
to hold at finite amplitudes. For small amplitude disturbances, the differ- 
ence between ( )* and ( )' becomes asymptotically small so that j- g F v c* 
— u'v' cos 2 (p is the usual horizontal momentum flux by the eddies. 

The Eliassen-Palm flux divergence defined in (2.9a) appears to be the 
most general form possible [It is essentially the same as in Tung (1982), sec- 
tion 5, but with several typographical errors corrected). Like in Andrews 
(1983), it is valid without restriction to small amplitude perturbations and 
quasi-geostrophic scaling. However, unlike in Andrews (1983), where adia- 
batic mean and eddy flows are assumed a priori, the present definition is valid 
for a general diabatic atmosphere. Furthermore, due to his use of the zonally 
averaged zonal momentum instead of the density- weighted zonal average, the 
eddy flux terms in his zonal momentum equation cannot be expressed in a 
pure divergence form as is done in our Eq. (2.8a), except under nonaccel- 
eration conditions. The present form permits the interpretation that eddy 
fluxes act to redistribute mean angular momentum (without net creation), 
and that when integrated over a volume bounded by a surface with no net 
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outward eddy fluxes, eddy forcing of mean flow vanishes. 

Nevertheless, the not-in-divergence form of eddy forcing terms of the 
zonally averaged momentum (which cannot be called an Eliassen-Palm flux 
divergence) is useful in a different aspect. It turns out that an exact relation- 
ship exists between this “Eliassen-Palm flux pseudo-divergence” , as we shall 
call it, and the flux of Ertel’s potential vorticity along isentropic surfaces. 
Since the flux of Ertel’s potential vorticity is not in general in a divergence 
form, no exact relationship exists between it and the Eliassen-Palm flux di- 
vergence defined in (2.9a). 

b. Eliassen-Palm pseudo-divergence. 

We shall now generalize Andrews’ (1983) expression for the eddy forcing 
term for the zonally averaged zonal momentum to include diabatic terms. 
Using the zonal average of u instead of the density-weighted average of u, 
we have instead of (2.8a) the following 

d d d 

— (p e ucos<p) + —(Vucos<p) + — (W u cos (p) - fV = 7, (2.8 b) 

where • 7 represents the eddy forcing of the zonally averaged absolute 

r$ 

angular momentum, and is given by 

1 n r - ~ i it 9 - \ i p ' e 

— □ • 7 = — tr cos <p—u' cos <p — If — — u cos <p)v' cos <p — 

Po oy dy pg 

cos<p + — p'gb'—ii cos <p. (2.9 b) 

dv Pg do 

By comparing (2.8b) with (2.8a), and noting that 

U — U = p'gU 1 cos<p /pg, 
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we see that 


□ • 7 = V • F - — [Vp'gu' cos(p/ p e ) 


d 


- ^(^PW'cos^/pfl) ~ — p' s u' cos <p, 


( 2 . 10 ) 


so that in component form, the “pseudo-divergence” is given by 


_ _ d d d 

~df v + d6 Je + di 7u 


where 


1 


J y = F v — — Vp' s u'cos <p = — Wcos <p 
Pi 


X 1 d 

7g — Fe — — Wp'gU 1 cos (p = — W'u'cosv? H — p*— — 
p e g aoX 


and 


J t = — Pau'cosy? 


( 2 . 11 ) 


The above expressions reduce to those of Andrews (1983) for adiabatic waves. 
The presence of the last time-derivative term in (2.10) is what makes Q • J 
not a pure divergence. [The same situation appears also in Andrews’ (1983) 
expression for adiabatic atmospheres]. In pressure coordinates, there is no 
distinction between (2.9a) and (2.9b), because the “density”, p p , in pressure 
coordinates is — and so there is no density perturbation, i.e. p' p = 0. In the 
present case, Q • 7 can be put into a divergence form only if it is integrated 
with respect to time over a lifecycle of a conservative disturbance. 

Conservation of absolute angular momentum: 
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A statement of the conservation of angular momentum appears to be 
more obvious if Eq. (2.8a) is combined with Eq. (2.6) to yield 

+ (2120) 
where La is the absolute angular momentum, with L defined as 

L = (u + fiacos <p\ cos 
When L is used instead of L, we have 

*! L+ *£ 1+ "» 1 “ D ' y (2 - 126) 

In the following we call complete the listing of the zonally-averaged equa- 
tions for use in our later computations. Some of the additional approxima- 
tions that will be introduced do not affect our above discussion on Eliassen- 
Palm fluxes, which are based only on the zonal momentum equation without 
approximation (beyond that made in deriving the primitive equations). 

2.4. Meridional momentum equation 

It appears, as we will show in section 7, that for a fast rotating planet like 
the earth, the dominant terms in the mean meridional momentum equation 
is expressed in the following balance 

d - , 

/u = — cos <p— <£ (2-13) 

oy 

where 

? = (20 H — ) sin<p. 

acos<£> 
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and $ is the Montgomery stream function. Eq. (2.13) becomes the geostrophic 
relationship if f is approximated by /. 

Within the same degree of approximation as (2.13), u can be replaced 
by u if (2.8a) is used instead of (2.8b). 

2.5. Hydrostatic relationship 


CpT = 0 §o* (2 - 14) 

2.6. Balanced wind relationship 

Combining (2.13) and (2.14), we find a diagnostic relationship between 
u and T. 

( u \ 8 8 

2 f2 H sin <p9—u = cos (p~~c p T. (2.15) 

\ a cos <p J 89 dy 

Again, u can replace u in Eq. (2.15). 

When the mean relative angular frequency u/(acos <p) is neglected when 
compared with the planetary rotational frequency 0, (2.15) becomes the 
thermal wind relationship 


fe 


8 _ 8 - 
-u = cos P - Cp r 


(2.16) 


(2.16) is a good approximation to (2.15) for Earth’s atmosphere. 
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2.7. Equation of state 


Using the hydrostatic equation 



(2.17) 


and the relationship between temperature and pressure from the definition of 
potential temperation 9 = T ^ P , we can obtain a relationship between 
density and temperature as 


Poo d (T\ c '! r 

96 g d9 \6 ) 


(2.18) 


To facilitate the taking of zonal averages, temperature is decomposed into a 
radiative equilibrium state, T e (y,0), and deviations from it, AT, i.e. 


T = T e (y,9) + AT, 


and it is assumed that AT is small compared to T t (y,9). (2.18) then yields 


_ _ Poo d 

g de 



(2.19) 


The assumption that AT/T e (y,9) is small appears to be valid in the strato- 
sphere even for finite amplitude disturbances. 


For future reference, we list the expression for the perturbation density 



3. Nonacceleration theorem. 


The mean, angular momentum equation (2.12a) is in the form of the 
so-called “generalized Eliassen-Palm theorem” (see Andrews and McIntyre 
(1976, 1978 a,b), Killworth and McIntyre (1985), although we have not shown 
that our — ~- t L has the same quadratic property as their fjA, the rate of 
change of the pseudo momentum. This is however not important for our 
purpose of obtaining nonacceleration theorems). If we let 

D d „ d \d 

Dt dt ad<p dO 

be the substantial derivative following the density-weighted mean circulation 
(M) = then Eq. (2.12a) is in the form 

a 

P*~L-V -F = 0. (3.2) 

[To the right-hand side of Eq. (3.2) one may add a term corresponding to a 
dissipative force that we have not included in our momentum equation.) 

The advantage of using isentropic coordinates becomes clear when we 
consider the situation under radiative equilibrium conditions. In the absence 
of diabatic heating, the thermodynamics equation W = q/t implies 

W = 0. (3.3) 


The continuity equation 


d d-d - 

51 f > + T» v + ae w 
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together with a polar boundary condition 


V = 0 at y = a or — a. 


then suggests that 


V = 0 


(3.4) 


at equilibrium. That is, there is no mean meridional circulation at radia- 
tive equilibrium. The same conclusion can be drawn in pressure coordinates 
only if geostrophic approximation is made in the zonal momentum equation. 
In nongeostrophic form in pressure coordinates, the presence of the advec- 
tion of mean temperature by ageostrophic mean circulation complicates the 
relationship between diabatic heating and the mean meridional circulation. 
Because of the direct relationship that exists in isentropic coordinates be- 
tween the mean meridional circulation and the presence of nonconservative 
process (such as diabatic heating), Eq. (3.2) can be written as 

d . - - 

p 6 — L-V'F = D, (3.5) 

where D should vanish in the absence of nonconservative processes (please 
note that a local time derivative is used in (3.5)). Using this so-called gen- 
eralized Eliassen-Palm relationship, we find that a necessary condition for 
nonacceleration in a conservative atmosphere is the vanishing of Eliassen- 
Palm flux divergence : 

V • F = 0 (3.6) 


Since 

□ •/ = V-F 
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under “nonacceleration conditions” , 


□ •7 = 0 (3.7) 

is also a necessary condition for nonacceleration. 

(3.7) is similar to the condition derived by Andrews (1983) and appears to 
be the most general version of nonacceleration theorems presently available 
in Eulerian coordinates. It is valid for finite amplitude waves satisfying 
primitive equations on a sphere. 

It should be pointed out that (3.6) is only a necessary condition for 
nonacceleration. Vanishing of the Eliassen-Palm flux divergence is in general 
not sufficient for nonacceleration. Furthermore, although adiabaticity, q = 0, 
implies V • F = 0 at steady state, as mentioned above, the converse is not 
necessarily true. That is, the vanishing of the Eliassen-Palm flux divergence 
does not guarantee that at steady state the atmosphere should be in local 
radiative (or radiative-convective) equilibrium, with q = 0. This observation 
has important implications concerning the presumed role of large-scale waves * 
in driving the atmosphere out of a radiative equilibrium. We will come back 
to this point in a later section. 
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4. Relationship between Eliassen-Palm flux pseudo 


divergence and the flux of Ertel’s potential 
vorticity. 


The relationship between the quasi- geostrophic Eliassen-Palm flux di- 
vergence in pressure coordinates and the flux of quasi- geostrophic potential 
vorticity is well known (Charney and Stern, 1962; Dickinson, 1969; Green, 
1970; Edmon et al, 1980). If the quasi-geostrophic potential vorticity can 
be assumed to be conserved along isobaric surfaces then a parameterization 
of the flux, and hence of the Eliassen-Palm flux divergence can be made 
in terms of the mean gradient of potential vorticity along isobaric surfaces. 
Such a procedure appears to be problematic because (a) the geostrophic as- 
sumption is too restrictive as it cannot account for mixing by gravity waves 
in the stratosphere, (b) quasi-geostrophic potential vorticity is in general a 
less-well conserved quantity along isobaric surfaces in the stratosphere than 
Ertel’s potential vorticity along isentropes, and (c) the quasi-geostrophic 
scaling requires that the cross-isobaric diffusion terms (such as K yz and K zz ) 
be neglected; this is perhaps not justifiable (see Tung (1984)) in pressure 
coordinates, except in limited regions in the mid stratosphere (see Newman 
et al (1985), personal communication). 

For these reasons it appears to be desirable to obtain a relationship 
between the nongeostrophic expression we have for the Eliassen-Palm flux 
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divergence and the flux of Ertel’s potential vorticity in isentropic coordinates. 
Then making use of the quasi-conservative nature of the isentropic potential 
vorticity, a parameterization of the eddy momentum forcing term (Eliassen- 
Palm flux divergence) can possibly be made. 


4.1. Ertel’s potential vorticity 


Ertel’s potential vorticity is defined as (see e.g. Pedlosky (1979)), 


n = 


S + f 

Pe 


(4.1) 


where 


f = 


d d 


— u - — ucos 
ox ay 


<P 


is the relative vorticity, and 


dz 1 dp 
Ps=P dd ~ ~gd9 


is the “density” in isentropic coordinates. [Recall also, y = a sirup, ^ = 


l a 

ocoa ip dX ' 


It can easily be shown, using (2.4), that the following identity holds for 
the perturbation (from its density-weighted average ) potential vorticity 

Q 

Pe n* =?'-(/- — u cosip) p' 9 /pj. (4.2) 

oy 

[Note that although (4.2) is of the same form as the linearized expression for 
p e IT', no small amplitude approximation is made in deriving it]. Using the 
identity that 

a'bt = a'f /, 
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we find that the northward flux of Ertel’s potential vorticity is given by 


pg(pgll*)v* cos<p = pep B R*v' cos <p 


= — p s v' cos>p — u'costp — (/ — — ucos<p)p' 6 t/'cos<p (4.3) 

The expression in (4.3) is to be compared with the Eliassen-Pakn flux 
pseudo-divergence, which is, from (2.9b) 


i—i -r _ . u 

M • / = —pev 1 cos <p— u' cos (p 
dy 


— if — — U COS ip) p'gV 1 COS ip 

dy 

' Q 

cosip + p' 9 0'~u cos <p 


(4.4) 


4.2. Adiabatic waves 

For adiabatic waves, (4.3) and (4.4) are the same, so 

□ • J — Pe{pe II*)v* cos <p (4.5) 

Eq. (4.5) provides an alternate diagnostic expression for the Eliassen-Palm 
flux pseudo-divergence (the net eddy forcing of the mean flow) which should 
be valid for time scales less than radiative damping times.. Since no small 
amplitude assumption has been made, (4.5) should probably hold even for 
large-amplitude “breaking” planetary waves observed in the stratosphere 
(McIntyre and Palmer, 1983, 1984; Leovy et al, 1985), to the extent that 
such waves tend to mix predominantly along isentropes. However, the ex- 
pression should not be expected to hold in the mesosphere in the presence 
of significant breaking gravity waves, which can mix momentum across isen- 
tropes. 
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4.3. Diabatic waves 


To give an estimate of the neglected terms in (4.5), we first write down 
the full expression without approximation, 


• d 


t: d 


□ • / = peipgU^v* cosip - pod' — u' cos tp + p' s 6'—u cos<p. (4.6) 

The last two terms in (4.6), denoted by A, will be estimated as follows. Since 

_ 2 ( 1 d V _ d , , d _ 

p *{7eM u ) " p °do u 

we have 


• d d ,• ( 1 d V 

A = -ped'—u 1 cos <p + p'gO'—ucosip ~ -p;6' ( ~ I 

The eddy cross-isentropic velocity O' is estimated using the Newtonian cooling 
expression (see Tung (1984)) 

•, <*0T' t 

6 = , where a ~ 5 days 

and the vertical derivative of u is estimated using the thermal wind relation- 
ship. Thus 


“:C„4-V 


Pe f 0 T P ad<p 2 
Using a value of T' ~ 10°C and a horizontal scale of a, we find 

^ ~ 10~ 6 m/s 2 
Pe 

This is about two orders of magnitude smaller than the Eliassen-Palm flux 
divergence, which is of the order of 


0 • 7 /pe ~ 10 *m/s 2 
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during winter, in the extratropics (Clough et al, 1985). 

Thus it appears that the relationship (4.5), which was found for adia- 
batic waves, should continue to hold for adiabatic waves to a good degree of 
accuracy. Thus, we have 

□ • 7 ~ p e {p e Tl*) v*cos<p (4.7) 

This is the diagnostic relationship that we have sought as a generalized ver- 
sion of the geostrophic relation commonly used. 
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5. A parameterization of Eliassen-Palm flux 
pseudo-divergence 

To incorporate the effect of diabatic as well as adiabatic eddies, we start 
with the nonconservative form of Ertel’s potential vorticity equation (see 
Pedlosky (1979)): 

= 5, (5.1) 

where 

is the diabatic source or sink of potential vorticity. When written in the form 
of (5.1), the transport of II obeys the same equation as that for (nonconser- 
vative) tracers. It turns out that because of the definition (4.1), the source 
term is proportional to the Ertel’s potential vorticity itself, i.e. 

S = (5.2) 


To obtain a parameterization of the flux of Ertel’s potential vorticity, 
we take guidance from the small amplitude theory to derive the form of 
functional dependence on the mean quantity. The perturbation form of (5.1) 
is 


Mh+ n *+ p'stM & "•+ M&ir 

(a) (b) (c) (d) (e) 


(f) 


(g) 


00 


0 ) 


(5.3) 
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As usual, the advection of TI* by the mean meridional circulation (t), 6) (terms 
(d) and (e)) is to be neglected when compared with the advection by the mean 
zonal wind (term (b)). Under the same approximation, term (h) should also 
be neglected as it is of the same order as or smaller than, term (e). An 
additional assumption is that the time scales of eddies are much shorter than 
that of the mean quantities, so term (c) should be neglected when compared 
to term (a), resulting in: 


d d » » d 


(5.4) 


where we have defined the displacements rj* and <t>* in the horizontal and 
vertical directions as 


<s + 


(5.5) 


Thus, 

(p fl IT)v* cos<p> = —v*(peTi*) cos 2 (p— ft 

°y 

d d 

- u*cos<p(p^*)^n + v*cos£>(p s ^<£*)n. 

(5.6) 

Let 



- 

d X 

p s K vv = u*(p S T7*) = Pe-^-r)' 2 

(5.7) 


PbK vS = v*(p 8 <f>*) 

(5.8) 

and 

A 



Q 

p e v E = v*{p e —<f>*). (5.9) 
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(5.6) then becomes 


a - 


a - 


{pe IT)t>* cosy? = -p e K vv cos 2 y?— f[ - p 6 K y e cosy?— IT 

ay do 


+ p$ve cos y?n 


(5.10) 


Thus (4.6) implies that 


1 d d 

— □ • ? = ~PeK yy cos 2 y?— n - p e K y$ cos y?— — IT + 
Pe dy v d$ 


(5.11) 


Since [] • 7 appears as the eddy forcing for the zonal mean absolute angular 
momentum L (see Eq. (2.12b)), (5.11) should also be rewritten in terms L. 
This can be accomplished by noting the following identity 


Pe n — / + ? = — ^-L 
dy 


(5.12) 


[In particular, the horizontal gradient of the Ertel’s potential vorticity is 
given by 


_ - _ _ _ a i a _ a 2 _ a 

+ —r--r-L 




dy p 9 dy 


a - 
dyP& 

ay 2 ~ ' Pe dy 


2fi a 2 % .(/ - ij' icos( p) T a 


= — -^(ucosy?)-/ 
a dy 2 


cos <pgp e dtnd 


Pc d _ 
ft din0 U 


(5.13) 


where 


Pc{0) = Poo fy 


cp/R 


Using the parameterization in (5.11), the mean zonal momentum equa- 
tion becomes 


(5.14) 


' dype dy 


- K yv cos 2 ippe-^-^-^-L + K v0 cos<pp 6 -^^--^-L. 


d_}_d_^ 

' dd p e dy J 


The coefficients K vv and K v0 for the diffusion of mean angular momentum 
turn out to be the same as those for the diffusion of conservative tracers as 
described in, for example, Tung (1982; 1984) and Ko et al (1985). Thus it 
appears that an assessment of the mean angular momentum budget or of 
the Eliassen-Palm flux divergence could also possibly afford us a preliminary 
look at the magnitude and latitudinal distribution of the diffusion coefficients 
used in zonally averaged model of tracer transports. These quantities are 
difficult to obtain by direct evaluation from data on transient waves in the 
atmosphere. Before we give a preliminary estimate of K vv in Part II, a few 
remarks concerning Eq. (5.14) are given in the following. 

(i) It is well-known that angular momentum L is not conservative and 
so cannot be treated in the same way as a conservative tracer. If L 
were a conservative quantity, i.e. if it were to satisfy 

d 


dt 


L = 0, 


instead of the actual equation 


1l = — — 

dt ad X 




which we have used, then the procedure outlined in Tung (1982) for 
conservative tracers would have given, for example, for the first term 
on the right-hand side of Eq. (5.14), 


cos 2 <Pi- L 


d . 
dy 
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instead of 


„ 2 - u A U T 

K vv cos <ppg — ~-^-L 

oyps dy 

actually appearing in Eq. (5.14). The difference should be attributed 
to the nonconservative nature of L, and not to the “quasi-adiabatic 
eddy approximation” used to obtain (4.6). 

(ii) If we assume that the vertical to horizontal gradient of the mean Er'tel 
potential vorticity is of the order of the ratio of a scale height to the 
radius of the earth, then estimates given in Tung (1984) would suggest 
that the magnitude of the K yi term in Eq. (5.14) (or (5.11)) is about 
10% of the magnitude of the K yv term in that same equation in the 
winter stratosphere in isentropic coordinates [These two terms are 
comparable in pressure coordinates under the same circumstances]. 
Therefore the right-hand side of Eq. (5.14) is dominated by the K vv 
term, during winter, 


Jt l + ( { + cos + e Te l = (5 - 15) 


(iii) The form of the right-hand side of Eq. (5.14) suggests that large-scale 
eddies act to (more or less) diffuse mean absolute angular momen- 
tum. It is perhaps important to note here the difference between the 
diffusion of mean absolute angular momentum and the diffusion of 
mean zonal velocity. The latter mechanism has sometimes been used 
as a parameterization of the eddy forcing term in the mean zonal mo- 
mentum equation. It can be shown that the right-hand side of (5.15) 
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can be written as 


K, 


d i d 9 


yy COS <PP$ _ ^ L 

dy Pe oy 


KyyC OS 2 <P 


d 2 

— (uc° s <p) + 


/(/- £(“ cos v?)) r d ( Pe d 


cos <pgp e dtnO \T e d£nO 


u - 


2Q 

a 


(5.16) 


If the mechanism of zonal velocity diffusion were adopted, the last 
term in (5.16), 

, 2fl 

— K vv cos <p — (5.17) 

a 

would have been absent. The presence of this term, which is propor- 
tional to the planetary rotational frequency, has important physical 
implications. This term is negative independent of the sign of mean 
gradients provided that the flux of Ertel’s potential vorticity is down 
gradient (so that K vu > 0). Therefore it provides a source of easterly 
relative momentum. The first two terms in the right-hand side of 
Eq. (5.26) act to diffuse relative angular momentum, thus serving 
to smooth horizontal and vertical gradients of the zonal wind u, but 
providing no net source of mean relative momentum. The presence 
of (5.17) serves to decelerate the westerly jet, balancing the Coriolis 
acceleration associated with the poleward mean meridional flow that 
exists in the stratosphere in the winter hemisphere. Without the term 
(5.17), the stratospheric westerly jet would reach unrealistically large 
velocities. By thermal wind, the temperature near the winter pole 
would be too low. 

The simple parameterization of Rayleigh friction sometimes used (Leovy, 
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1964; Holton and Wehrbein, 1980): 

~U-^ = -k r u, 

Pe 

also provides the needed easterly momentum source in the westerly 
region. In this sense, Rayleigh friction is probably a better parame- 
terization than the one based on diffusion of u, i.e. 

• / ~ fiV 2 u. 

Pe 

(iv) If the same balance of terms as in Eq. (5.15) were to prevail in the 
easterly region in the summer hemisphere, the presence of the easterly 
momentum source (5.17) would have led to an unrealistically strong 
easterly jet. However, the presence of the same easterlies prevents 
most of the stationary planetary waves from penetrating above 30 km, 
so that in this region, K vv should be very small. Also the estimate of 
the ratio of K v $ to K vv terms mentioned earlier for winter stratosphere 
probably does not apply to the summer hemisphere. If the eddy 
forcing term is dominated by the effect of diabatic gravity waves with 
short vertical scales, the K v g term in Eq. (5.14) may even dominate 
over the K vv term in the easterly region. Since K v g can have either 
sign, the problem of easterly acceleration in a easterly region does 
not necessarily arise. 

(v) Strictly speaking, one should not treat (5.11) or (5.14) as a parame- 
terization of the eddy forcing of the zonal mean flow unless the depen- 
dence of the coefficients K vv and K v $ on u or L is also known. This 
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situation for the momentum equation should be contrasted with that 
for the tracer transport equation, where the K ' s are independent of 
the. concentration of the minor species being transported. Neverthe- 
less, Eq. (5.11) or (5.14) is useful for diagnostic purposes as one can 
in principle deduce from the momentum budget the same set of K ' s 
that are used in tracer transport equations. 
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6. Relation between Eliassen-Palm flux diver- 


gence and diabatic heating. 

The often-cited linear relationship (see WMO/NASA (1985), Chapter 7) 
between Eliassen-Palm flux term and the diabatic heating is obtained only 
when the geostrophic approximation is applied to the mean zonal momentum 
equation, which becomes, 3 


-fv= 


( 6 . 1 ) 


The steady state continuity equation is, from (2.6), 


r/ + ie*= a - 


( 6 . 2 ) 


[Incidentally, it has been shown in Tung (1982) that (6.2) is a good approx- 
imation to (2.6) even without the steady-state assumption]. The thermody- 
namics equation is, in isentropic coordinates, from (2.7): 


W = q/T 


(6.3) 


Given the eddy forcing, □ • J, Eq. (6.1) determines the meridional velocity 
V locally at every point where the eddy forcing is specified. The vertical 
velocity W is calculated from V through the continuity equation (6.2), given 
a boundary condition on W at some level. The thermodynamics equation 
then relates W to the diabatic heating field q. 

3 Note that under the geostrophic approximation, V F and Q 7 can be used interchangeable 
in Eq. (6.1) 
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The physical implications of these relations have often been taken to 
be that (a) it is the eddies that are responsible for driving the atmosphere 
away from local radiative equilibrium, and (b) the degree of deviation from 
radiative equilibrium is directly proportional to the strength of the Eliassen- 
Palm flux divergence, the last point being inferred from the relationship 

|(3 / r ) = | ; (i D .r), (6.4) 

obtained by combining (6.1), (6.2) and (6.3). 

There appear to be two points of caution that we wish to inject con- 
cerning the interpretation just mentioned. First, as mentioned above, W 
and hence q , are influenced by boundary conditions as well as in situ eddy 
forcing, and so the above results apply only in the absence of significant 
boundary forcing in W. Second, geostrophy is not a good approximation 
when applied to the zonally averaged momentum equation in the east-west 
direction (It is a better approximation when applied to the meridional mo- 
mentum equation, yielding the thermal wind relationship). Results on zonal 
mean circulations using the primitive equations are known in some cases to 
differ significantly from those obtained using geostrophic equations. 

This model difference can be understood by examining the relative im- 
portance of various terms in the zonal momentum equation: 
d d 

V — (u cos<p) + w— (u cos <p) - fv =□• 7 (6.5) 

ay ov 

In geostrophic approximation, the nonlinear advection terms, the first two 
terms in Eq. (6.5), are assumed to be much smaller than the Coriolis term, 
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the third term in Eq. (6.5). Consequently the balance is between the Coriolis 
term and the eddy forcing term on the right-hand side (see Eq. (6.1)). As 
□ • 7 is made smaller and smaller, the Coriolis term necessarily decreases 
in importance according to the geostrophic balance (6.1). This eventually 
leads to the breakdown of the geostrophic approximation as the neglected 
relative acceleration terms are no longer small when compared to the Coriolis 
force. The unanswered questions are: What is the threshold magnitude 
for eddy forcing below which the geostrophic balance (6.1) no longer holds, 
and the “nonlinear nearly inviscid regime” (Held and Hou, 1980) applies? 
Which regime is our atmosphere in? An assessment of the magnitude of 
[] • J and the validity of geostrophic approximation as applied to the mean 
zonal momentum equation are necessary for an understanding of these issues. 
The diagnostic assessments will be given in Part II. In this section, some 
conceptual issues concerning the state of the atmosphere in the absence of 
eddy forcing will first be addressed. 

The steady state zonal momentum equation can be written as 

J[*,L]=X, (6.6) 

where is the streamfunction of the meridional circulation defined by 



The Jacobian operator appearing in (6.6) is defined by 

J [A, B] = A v Be — A$B V . 

In (6.6) X represents the large-scale eddy forcing plus small-scale dissipative 
terms, i.e. 

AT = D'J + V-mVZ (6.8) 

where is the coefficient of small-scale mixing or molecular diffusion. In 
the absence of zonal asymmetries, X is represented by only the (molecular) 
diffusion. 

6.1. Inviscid symmetric states. 

If the hypothetical atmosphere is symmetric ([] • J = 0) and inviscid 
(fj. = 0), (6.6) becomes 

j[$,I,]=0, (6.9) 

which implies that lines of constant absolute momentum should coincide with 
the streamline. This is simply a statement that the inviscid symmetric cir- 
culation is angular- momentum conserving. There is no constraint from the 
momentum equation that the solution should be in local radiative equilib- 
rium: 

$ = 0 ( 6 . 10 ) 

(which implies, from thermodynamics equation, 

9 = 0). 
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This situation should be seen in contrast with the geostrophic version of the 
zonal momentum equation (6.1), which in the symmetric inviscid case is 

fV = 0, 

which leads to radiative equilibrium q = 0 in the absence of boundary forcing. 
The exact solution to (6.9) is 

* = G(L)- (6.11) 

the functional form of G is to be determined by boundary conditions. If the 
horizontal boundaries are in radiative equilibrium, i.e. 

W = 0 at 0 O and 0i, (6.12) 

we have, without loss of generality 

$ = 0 at 0 O and 0 X . (6.13) 

If every line of constant L touches either the upper or lower boundary, 
then (6.11) becomes, when evaluated using (6.13), 

$ = G = 0 (6.14) 

This is the radiative equilibrium solution, with q = 0. The zonal momentum 
is related through the balanced wind (or thermal wind) relationship to the 
radiative equilibrium temperature T e satisfying q(T e ) = 0. 

However, if in the interior of the domain there exists a region of closed 
contours of L , then there will exists “modon”-like solutions for which G is not 
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necessarily zero. There is in general a meridional circulation in such a region 
(an out-of-radiative equilibrium situation). Also, inside the closed circula- 
tion region the zonal momentum is no longer constrained by the radiative 
equilibrium temperature, and can thus in principle attain any distribution. 

The inviscid symmetric circulation is nonunique because neither the 
shape of the closed region, the functional form of G, nor the distribution 
of L inside this region is determined without the imposition of additional 
constraints. A common misconception has been that the radiative equilib- 
rium solution (6.14) is the only solution in the absence of eddies or friction. 
[Ironically, even Held and Hou (1980), who have established that a realistic 
looking Hadley circulation can exist even in the absence of large-scale eddies, 
chose to interpret their deduced circulation as due to the presence of small 
but nonzero viscosity as a singular perturbation to the inviscid solution, 
which they assumed to be in radiative equilibrium.] 

However, that the inviscid steady state recirculating solution can be 
nonunique is well known in other areas of fluid mechanics (see e.g. Batche- 
lor (1956) in the case of vorticity distribution in steady laminar flow, Stern 
(1975) for “modons”, and Tung et al (1982) for internal waves of “permanent 
form”). There axe at least two ways of removing the arbitrariness. These 
depend on whether the steady inviscid solution is treated as the limit of 
fj, — ► 0 and then t — * oo, or as the limit of t — ► oo and then n — > ► 0. In the 
first approach, one considers the “quasi-steady” limit of an inviscid flow (see 
Benney and Ko (1978) and Tung et al (1982)) and G can take any arbitrary 
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form depending presumably on initial conditions (if these are stable). In the 
second approach, one considers the “nearly inviscid” limit of a steady flow 
(see Batchelor (1956)). The presence of even a small amount of viscosity 
(which however acts for an infinitely long time) can place severe constraints 
on the distribution of G inside a recirculating region. 

6.2. The nearly inviscid limit. 

The presence of a small amount of viscosity eliminates the radiative equi- 
librium solution (6.10). More precisely, (6.10), or an order n modification 
of (6.10), cannot remain as a solution to the viscous equation even in the 
limit as /z — ► 0. The argument was given by Held and Hou (1980) using a 
version of Hide (1969)’s theorem and assuming that the imposed radiative 
equilibrium temperature decreases poleward for all heights. 

When a small amount of viscosity is present, the inviscid recirculating 
solution mentioned above it remains valid to leading order except in certain 
“boundary layers” , in which it may be modified in different ways depending 
on the boundary condition. If the imposition of viscous boundary conditions, 
such as no-slip, is inconsistent with the inviscid solution, a (thin) boundary 
layer is needed. Inside the boundary layer, the effect of viscosity is important 
and (6.9) no longer holds. Here, the constraint from (6.11) that streamlines 
follow lines of constant absolute angular momentum is broken. Consequently, 
the streamlines can be closed, as demanded by mass conservation, while lines 
of constant angular momentum can be open and intersect the lower boundary, 
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as required for example by the no-slip condition, u = 0 (i.e. L = ftacos 2 <p). 

Outside the boundary layer, the interior solution is almost inviscid, i.e., 

$ ~ G(L) (6.15) 

except now G is to be determined by matching to the boundary layer solution, 
and not by evaluating (6.15) at the boundary itself (where ^ = 0). 

The detail matched asymptotic solution will be presented in a separate 
paper. Here we only wish to emphasize the fact that the interior nearly 
inviscid solution (6.15) can be very close to one of the inviscid recirculating 
solutions (6.11). And it is therefore more natural to interpret the nearly 
inviscid solution to be a small perturbation (at least in the interior ) from 
a recirculating inviscid solution, 4 instead of it being a singular perturbation 
(with order one change) from the radiative equilibrium solution as suggested 
by Held and Hou (1980). 

The conclusion that one can draw from the above discussion appears 

to be that eddies (large or small scales) are not “needed” to “drive the 

atmosphere out of radiative equilibrium”. Nevertheless, sufficient viscosity 

may be "needed” to maintain s ome form of stability for these nearly inviscid 
4 In this case the inviscid solution can be taken as a modon-like recirculating configuration 

with the lower interface brought arbitrarily close to the lower boundary. The lower interface 

will in general be a region of discontinuity in which the lines of constant L will complete the 

closed contour. The addition of a small amount of viscosity will introduce a thin boundary 

layer around the discontinuity and, depending on the lower boundary condition, lines of 

constant L do not necessarily have to close via this boundary layer. 
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steady symmetric circulations. 



7. A modified scaling appropriate for zonal 


mean circulations 

In sectiion 6, two extreme circulation regimes are discussed. One, the 
inviscid or nearly inviscid regime, pertains to the mean circulation in the 
presence of no or small eddy momentum forcing. The other, the geostrophic 
regime applies away from the equator when eddy momentum forcing is 
strong. It was clear from the discussion that the geostrophic approxima- 
tion, obtained using Rossby number scaling, is of only limited validity when 
applied to different zonally averaged circulation regimes. 

In this section, a somewhat different scaling based on the smallness of 
the ratio of mean meridional and zonal velocities is proposed. It turns out 
that the result is more general and can cover both extremes of circulation 
regimes mentioned above. 

Since the following discussions are not restricted to isentropic coordinates 
only, the more commonly used equations in pressure coordinates will be used 
first, with modifications for isentropic coordinates given later. 


7.1. Implications of Rossby number scaling 

Implicit in scalings based on the smallness of Rossby number 

( 7 . 1 ) 

where u 00 is a typical value of the horizontal velocity, 5 i a typical horizontal 
5 Note that Rossby number scaling does not distinguish between the magnitudes of merid- 


*o|<= “ oo/ * 


20 sin (p 0 ’ 
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length scale of the flow, and 2H sin p 0 a typical (midlatitude) value for the 
Coriolis parameter /, is the assumption that the gradient of relative momen- 
tum be much smaller than the gradient of planetary momentum, resulting in 
the approximation 

D d 

— • L ~ v— (tlacos 2 p) = -fv (7.2) 


where 


and 


D_ddd _ ( Poo. 
m = ai + v T y +w Tz' z = >- w = 



L = u cos p + 0 a cos 2 p. 

(7.2) is the so-called geostrophic approximation, and by its nature necessarily 
excludes the kind of absolute-angular-momentum-conserving flow cited pre- 
viously and by Held and Hou (1980), which, by approximately conserving 
absolute angular momentum, L, generates in its northward branch west- 
erly zonal flows as the planetary angular momentum, Ha cos 2 p, decreases 
northward. Therefore, a characteristic of such a flow is that the variation 
(with latitude) of the relative zonal angular moemntum is comparable to the 
variation in the planetary angular momentum (even though u/(flacos<£>) is 
small). In other words, the Rossby number for such a flow should be order 
one and so cannot be used as an asymptotic parameter for this case. 


In our present problem of calculating the zonal mean circulations, the 

scale of the mean circulation should be obtained as part of the solution 

instead of specified a priori. 
ional and zonal velocities. 
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7.2. An alternate scaling 


We propose an alternate scaling based on the assumption that for a 
rapidly ro'tating planet, such as the Earth, the zonally averaged flow should 
predominately be in the zonal direction. Therefore, the ratio 


_ u oo 
Uoo 


(7.3) 


should be small. In (7.3), t) 0 o and u 0 o are the typical magnitudes of the zonally 
averaged meridional and zonal velocities respectively. [Note that there is no 
assumption that | - | be small] 

If, as is usually done, one nondimensionalizes both horizontal velocities, 
u and u, by the same uoo, then the following scalings should result: 


u, = u/uoo = Uq + eui + e 2 u 2 + 


and 


V, = v/uqq = £Vi + £ 2 v 2 + 


(7.4) 


Let t y and t z be the typical meridional and vertical length scales, then 
based on the continuity equation, the mean vertical velocity w should have 
the scale (f^voo), and so the usual dimensionless quantity 


u>. 


= w 



/uoo 


should be of order e, i.e. 


w, = ewi + e 2 w 2 + • • • . 


(7.5) 
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For mean transports in the meridional plane, time t should be scaled by 
ty/vaoi a nd so we define, 




Starting with the zonally averaged zonal momentum equation, which we 
write in the following form 

n 

(7.6) 


7T L = 

Dt 


where X includes the Eliassen-Palm flux divergence (if in pressure coordi- 
nates, the residual quantities are used for v and w), and dividing by 
ul 0 /£ y , we have 


d , , a . d 

e aF + (£5 ‘ + ->^ + (£ *‘ + -fe 


((u 0 + eui H ) cos <p] 


where 


— (eiii + • • •)/* = eXi + e 2 X2 + 


y t = y/t v and /, = f{i y /u 00 ). 


(7.7) 


In (7.7), we have also defined and scaled 


X, = X/{u 2 0 /l v ) = eXi + e 2 X 2 + 


(7.8) 


[The scaling with X . being order one would have violated our assumption of 
| v/u | small. This exceptional case will not be discussed here]. Therefore, 
to leading order in e, (7.7) yields 


d d d 

df + + 1 ( fioCOSV?) " Ui/ * = 


dz t 


(7.9) 
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Note that both the relative acceleration and Coriolis force are retained in 
this approximation. In the ^-direction, we start with the following zonally 
averaged meridional momentum equation 


d _ - ia= 


(7.10) 


where Y is the eddy forcing for the j/-momentum and / = (2fl + a ^ Sip ) sin<p. 
The same procedure as before yields 


d . d . d 


[eui + • • • 


+ (/o + + • * •)( u ° + + • • •) 


= -cosp— ($o ) + (eTi + e 2 Y 2 ■+ ), 

&y* 


(7.11) 


where 


/o — /* + Uq 



tan <p 


fi = 



tan <p. 


To leading order in e, (7.11) is 


/ouo = - cos<p^-$ 0 , 

oy . 


(7.12a) 


which is the balanced wind relationship. Modification to (7.12) from the 
mean meridional advection terms arises at order e 2 . Thus, the balanced 
wind relations (7.12) holds to high accuracy provided the eddy meridional 
momentum forcing is not too strong (i.e. if Yi = 0). 
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Eq. (7.12) reduces to the geostrophic relation if one further assume that 


the Rossby number based on the radius of the earth a, 

D t _ «/cos(p 


(7.13) 


is small. This is generally true for a rotating planet whose planetary rota- 
tional frequency, fl, is larger than the relative angular frequency, ti/(acos <p), 
of the flow on it. [Note that this ratio does not involve the unknown length 
scale £ v .) Thus a special, relevant case of (7.12) is 


f t u 0 = - cos<p-— $ 0 
ay , 


(7.14a) 


This geostrophic relationship for the zonal flow is obtained without the Rossby 
number scaling, i.e. the smallness of the Rossby number defined in (7.1) is 
not assumed. 


Assuming = 0, the next order in (7.11) is 

(/, + 2u 0 [ — tan <£>) Ui = — cosv9-^-$i, (7.126) 

\a J dy t 

and if furthermore (7.13) is small, (7.12b) becomes again the geostrophic 
relation: 

Q 

f,u 1 = -cos<p-^— (7.146) 

ay . 


The energy equation, 



9 . 

T ’ 


when rewritten in log-pressure coordinates, becomes 



where Too is the reference temperature associated with the reference scale 
height, H 0 o = RToo/g, used in the definition of z. The zonal mean of the 
above equation is 

w, T+a hi =Q+ ^ (7 - 15) 

where Z contains the eddy heat fluxes. We shall assume that (t/,u>) is the 
residual mean circulation, so that Z is not of dominant importance in the 
energy equation. This assumption will be reflected in our scaling for Z to 
follow. 

We choose 

T* = T/Too = To + eTi + • • • 

Z * = Zi y /(u 0 oToo) = zZ\ + e 2 Z2 + • • • (7.16) 

and Q * == Q£v/{u 00 T 00 ) = Oo + &Qi + • • • 

For our present purpose, we shall assume that the mean radiative heating 
rate Q is a function of the mean temperature, viz, 

0 = Q(T). 


Therefore, 

Q. = Q.(T.) = Q.(%) + + ■ ■ • (7.17) 

ol 0 

and thus 

Oo = Q.(Tb) (7.18a) 

Q 

Qx = ajrQ*(T Q ) T u etc. (7.186) 
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Eq. (7.15), when divided by (iiooToo)/£y, becomes 


d . d d 

£— + (euj H ) — + [ewi + • • •) -r— 


dt< 


dy * 


dz. 


[Tq + eT\ + • • •] 


+( £U, i d )(To + zTi d ) 

-tOO Cp 


— Qo + sQ 1 + • • • + + 


(7.19) 


The leading term in (7.18) is 


Qo = 0, 


(7.20) 


which defines the leading order temperature field via (7.18a) as 


% = T e , 


(7.21a) 


where 


§,(T e ) =0, 


(7.216) 


In other words, To is the radiative equilibrium temperature T t satisfying 
(7.21b). At the next order in e, (7.19) yields 


' d d _ d \ - - 

tt— + Vi — h tZ>i~ — T 0 + w\Tq 

.dt, dy. dz„) 



= Oi + Zu (7-22) 

which together with (7.18b), should determine 7\, the mean temperature 
change (from radiative equilibrium) induced by the mean meridional circu- 
lation (vi,tDj). 
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7.3. Summary of scaled mean equations in isobaric co- 


ordinates 

In dimensional form, we write the zonal mean temperature and zonal 
flow as 

T = T e {y,z) + AT (7.23) 


u = u e (y, z) + Au 


(7.24) 


where T e is the local radiative equilibrium temperature determined from 


Q(Te) = 0 


(T.25) 


And u e is related to T e through the balanced wind relationship 

2(0h - — ) s\n<p-^-u — cos<p Q r~~f (7.26) 

\ a cos <p ) oz Hqq 

(Substitute u e for u and T e for T in (7.26)) (7.26) is obtained by combining 
the hydrostatic relation 

a 

(7.27) 


= RT t 

oz 


with the balanced equation 


f _ \ Q 

20 -) — — ) sin<p • u = — cos 92 — (7.28) 

^ a cos <p J oy 


which is correct up to and including order e. 


The zonal momentum equation is 


'd d d\ , ^ 

— + V— + w— j u e cos <p — fv = X 


(7.29) 
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Knowing u e from (7.25) and (7.26), this is a linear equation for the meridional 
circulation (u,tD) assuming that the eddy momentum forcing, X, is known. 
v and w are related to each other by the continuity equation, which is 

d \ d 

— (0 cos ip) + — (p 0 ( 2 )to) = 0, (7.30) 

oy p 0 z oz 

where Pa{z) = po(0)e -z ^ 00 . 

Knowing v and w, the energy equation 

(£i + % + T ' + 1 

+ ■ AT (7.31) 

then yields the dynamically induced part of the mean temperature A?. And 
from A?, the balanced wind relation gives diagnostically the dynamical cor- 
rection to u e . Note that the solution procedure outlined above for the scaled 
equations is linear in every step, while the original set of mean equations is 
nonlinear. 

Note also that the above set of equations includes the geostrophic equa- 
tions as a subset; one can make the additional assumption of small Rossby 
numbers, Ro |<, and the geostrophic equations will result from the above set. 
It also includes the nearly inviscid regime mentioned in section 6; even for 
small or zero eddy forcing X, a mean meridional circulation (u,u>) and an 
out of equilibrium temperature distribution T e + A T can still be the solution 
of our scaled set of equations, given appropriate boundary conditions. 
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7.4. The scaled mean equations in isentropic coordi- 


nates 

The set of scaled mean equations in isentropic coordinates is the same 
as stated in section 2, except 

(i) with the radiative equilibrium temperature T e (y, 6) replacing T in the 
static stability parameter appearing in the thermodynamics equation 
( 2 . 7 ). 

(ii) with u e replacing u in zonal the momentum equation (2.8a) and u in 
(2.8b), and similarly in (2.12a) and (2.12b). 

(iii) with q(T) approximated by q = ^ q(T t ) • ( T — T e ). and 

(iv) under the present scaling, there is no difference between V • F and 
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8. Conclusions 


A formulation is given for a nongeostrophic theory of zonally averaged 
circulation. Motivations for extending the commonly used geostrophic ver- 
sion of the general theory are many; some of them are 

(a) Geostrophy does not hold at and near the equator and so is not 
uniformly valid in a model of zonal mean global circulation. 

(b) Although there is some consensus that the mean zonal velocity is 
in geostrophic balance, there is no a priori justification for applying 
the geostrophic approximation for the zonally averaged north-south 
velocity. This is the case even if the nonaveraged north-south veloc- 
ity can be considered as in geostrophic balance with the east-west 
gradient of pressure (or height) field. 

(c) In the eddy forcing terms for themean circulation equatorial wave 
and gravity waves, which are important components of the forcing, 
are filtered out in the geostrophic version of the theory. 

(d) The constraint of geostrophy is so strong that one is prevented from 
addressing some fundamental problems concerning zonally averaged 
circulations using the set of geostrophic equations. For example, it is 
probably not appropriate in using the set of geostrophic equations to 
deduce the state of the atmosphere in the absence of zonal asymmetry, 
because while these equations would lead one to conclude that the 
symmetric atmosphere is extremely close to radiative equilibrium, 
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while the original primitive equations suggest otherwise. 


A major hindrance in the past in adopting the more general(i.e. non- 
geostrophic) formulation of zonally averaged circulations lies not in technical 
difficulties associated with the formulation but in clear interpretation of the 
role of various eddy forcing terms. This problem is largely circumvented 
when the zonal averaging is performed on isentropic surfaces, as it is done 
in the present work. The role of eddies in the forcing of the mean flow is 
much better defined in isentropic coordinates without the presence in isobaric 
coordinates of mean temperature advection by the ageostrophic circulation. 

The more general formulation suggests that a hypothetical atmosphere 
can be in different circulation regimes depending on the strength of the 
eddy Eliassen-Palm flux divergence. When the eddy strength is strong, the 
mean circulation may be in the geostrophic regime away from the equator. 
However, when the eddy strength is weak, the circulation may be close to 
the “nearly inviscid nonlinear regime” of Held and Hou (1980), which is not 
a continuation of the geostrophic regime as the eddy strength is reduced. 

In Part II, we will address the question of which regime our atmosphere 
is in. Furthermore, we will use our present isentropic formulation and the 
derived relation between Eliassen-Palm flux pseudo divergence and the flux 
of Ertel’s isentropic potential vorticity to deduce the isentropic diffusion co- 
efficient, K vv that is useful in tracer transport studies. 


57 



ACKNOWLEDGMENT 


Part of the research was carried out while the author was a Guggenheim 
Fellow on- sabbatical leave from M.I.T. The generosity of the John Simon 
Guggenheim Foundation is gratefully acknowledged. The research is spon- 
sored by National Aeronautics and Space Administration under grant no. 
NAGW-798. 


Appendix A. 

The zonally averaged momentum equation in isentropic 
coordinates on a sphere. 

Starting from the zonal momentum equation: 


8 u du v 8 ■ 8 

— u H 1 u + 9 — u 

dt a cos pd\ a dp 88 


- 20 + 


u 


sm pv = 


a cos p J a cos pd\ 

and the equation for the conservation of mass: 


$ 


(A.l) 


8 Id 18. . 8 . •. . . 

— ( Pev cos p) + — (p*0) = 0, (A.2) 


dt a cos pd\ 


a cos <p dip 


one has 


8 Id Id 

— p fi u+ p«uuH — — (p^uvcos p) + — (potifl) 

dt a cos pd A a cos 2 p dp 88 


Pe d 

20 sin ppev = — $ 

acosp d A 


(-4.3) 
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Taking the zonal average: ( ) = / 0 2r ( )d \ , one obtains the zonally averaged 
momentum equation: 


Q 1 Q Q 

—pju + — — (pjuv cos 2 p) + —( p e u'e ) - 2fi sin p~p$v 

at a cos 2 p dp ad 


Pe d 


a cos p d\ 


$ 


(A.4) 


We let, as in Tung (1982): 


U = p$u,V = pgvcosp, and W = 


(A.5) 


be the zonally averaged mass circulation, and define, as in Gallimore and 
Johnson (1981a) the density-weighted zonal average 

u = Wu/pg = U/Pe, 
and deviation from zonal average 


u = u — u. 


(A- 6) 


Then 


pguv cos 2 p = pi [u + u* 


V 

— + V* COS ip 

Pe 


COS ip 


= uV COS p + PgXL*V * cos 2 p + ~Pe u * COS pV + XL COS 2 PPeV* 

Pe 

The last two terms are zero from definition (A. 6). Eq. (A.4) now becomes 


—p e xL -i — — — 

at a cos ‘pap 


Q 

(u cos pV’') + — (uW^j — 2flsiny?V /cos p 


$ ) Pe ~ 7 _ — {pgXL'v* cos 2 p) - — (pinTfl*) (A. 7) 


acosp dX 


a cos 2 p dp 
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The right hand side of (A. 7) can be written as 


1 - - 1 
-V • F = 


cos p cos p 


1 d „ d „ 

■Fy + —Fg 


a cosp dp dd 


where 


F y = —p e u*v* cos* p, 


and 

— r 1 5 ' 

Peu *9* p —$ . 

g a cos pd A 

In order to write the first term on the right-hand side of (A. 7) into a vertical 
derivative term, we have made use of the hydrostatic equation 

d 

gPe ~ dO P 




and the relationship between the Montgomery streamfunction $ and pressure 



Finally, by multiplying both sides of Eq. (A. 7) by cosp, we obtain 


d \ d d - 

— [p e u cos p) H — (u cos pV ) + — (u cos pW ) — 2fl sin pV = V • F 


dt 


a cos p dp 


(A- 9) 


Next, the mean momentum equation expressed in terms u (instead of u 
as in Eq. (A. 9) is derived. We start with Eq. (A.l) and perform the zonal 


average, 


d v d d 

— (ucos p) H — — — ucos p + 9— u cosp — 2Cls'mpv cosp = 0 (A. 10) 

dt a dp do 
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We let v = v + u* , but u = u + u', then 


v d v d v* d v* d 

— — - u cos tp = - — u cosip 4 — — u' cosp H — — ucos p 

a dp a dip a dp a dp 


Since 


v* = v' - p'dv’/pb 
from the definition of v*, and so 


a*bf = a'bf , 


we have 


v d 


v d 


1 


-1 d 


- — — u cos p = - ——u cos p H ——u'cosp — — pit/-— -u cos p. 

a dp a dp a dp p e a dp 


Similarly, 


- d • d • d 1 — ~ d 

9-^zU cos V ~ cos V + cos V ~ ~P'» d '^2 u cos 'fi- 

at) at) at) po at) 

Substituting 


v = v - p'gv'/pe 


into (A. 10), we find 


d . . v d J 9 _ _ . 

— (u cos p) H — — — cos p + 0— u cos p — 2\l sm<pv cos p 
dt a dp d6 


v' d • d 

= — — u' cos p — 9'— u'cosp 

a dp dd 


— 2 fl sin p — 


r d 


a cos p dp 


ucosp p'gv'cosp/pe + — p'gO'— ucosp (A.ll) 
I p B dd 


Using the zonally averaged continuity equation (from (A. 2)): 


d - 
Ft ! " + 


1 d . . . d 

a - \ psV cos V>) + 37 

a cos p dp do 



= 0 


(A.12) 
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we rewrite (A. 11) finally into the form 


where 


d . 1 d . _ . 

— (Pffii cos <p ) H — — ( p e v cos <pu cos <p) 

at a cos <p dip 

+ — (j>$6ucos<pj — 2fl sin ipv cos p = 0 - 7, 


1 u' 5 • d 

— □ • 7 — u'cos p — O'—u 1 cos <p 

Pe a dip do 


^ \ ^ ^ 

2H sin v? — u cos p'^' cos<p/p e + -p' 9 6'-—u 

a cosip dip J p e d8 


COi 


(A. 13) 


p (A. 14) 
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